Matter instabilities in general Gauss-Bonnet gravity 
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We study the evolution of cosmological perturbations in f(Q) gravity, where the Lagrangian is 
the sum of a Ricci scalar R and an arbitrary function / in terms of a Gauss-Bonnet term Q. We 
derive the equations for perturbations assuming matter to be described by a perfect fluid with a 
constant equation of state w. We show that density perturbations in perfect fluids exhibit negative 
instabilities during both the radiation and the matter domination, irrespective of the form of f(Q). 
This growth of perturbations gets stronger on smaller scales, which is difficult to be compatible with 
the observed galaxy spectrum unless the deviation from General Relativity is very small. Thus f{Q) 
cosmological models are effectively ruled out from this Ultra- Violet instability, even though they 
can be compatible with the late-time cosmic acceleration and local gravity constraints. 



I. INTRODUCTION 

Independent observational evidence for dark energy 
has stimulated the idea that General Relativity (GR) 
may be modified on large distances to give rise to a late- 
time cosmic acceleration A simple dark energy sce- 
nario constructed in this vein is so-called f(R) gravity in 
which / is a function of the Ricci scalar R |2j . Although 
there are some restrictions to the functional form of f{R) 
to satisfy both cosmological and local gravity constraints, 
it is possible to design viable models |3( that can be dis- 
tinguished from GR at least in the metric formalism of 
f(R) gravity. 

The f(R) gravity in the metric formalism corresponds 
to the so-called Brans-Dicke theory with a parameter 
^bd = in the presence of a potential of gravitational 
origin Q. One can generalize this to scalar-tensor the- 
ories with an arbitrary Brans-Dicke parameter ojbd- In 
fact it is possible to construct scalar-field potentials that 
can be responsible for the cosmic acceleration, while at 
the same time satisfying local gravity constraints Q. 
These models, including f(R) gravity, exhibit several in- 
teresting observational signatures such as the phantom 
equation of state @ , the modified matter power spec- 
trum Q , and the modified weak lensing spectrum (8[ . 

The Ricci scalar R is not the only scalar quantity which 
is used to change gravity, as we can easily construct other 
scalar quantities such as R^R^ and R ilvpa Rf il ' pa from 
the Ricci tensor R pl , and the Riemann tensor R^ vpcr @. 
However, for the Gauss-Bonnet (GB) curvature invariant 



g = R 2 - ARuuR^ 



1 ^pvpa 1 L 



(1) 



one can avoid the appearance of spurious spin-2 ghosts 
fiol fllj . If a scalar field qb with an exponential potential 
V((j>) = Vbe _A< ^ couples to the GB term [12], a scaling 
matter era can be followed by a late-time de Sitter so- 
lution for the exponential GB coupling F(4>) cx with 

n > a mini- 

While this GB coupling is well moti- 
vated by low-energy effective string theory (l5j | the joint 
likelihood analysis using observational data of big bang 



nucleosynthesis, large scale structure, and baryon acous- 
tic oscillations disfavors such a model, provided that it 
aims to account for dark energy [l6j . In addition the 
energy contribution coming from the GB term needs to 
be strongly suppressed for consistency with solar-system 
experiments pjl, [H|- This cannot be compatible with 
the requirement of cosmic acceleration today (see also 
Ref . [l9| ) j a * least in the presence of a kinetic term for the 
scalar field and some forms of the potential. The instabil- 
ity of tensor perturbations is also present in those models 
during the epoch of cosmic acceleration 0, [2(| Hl[ . 

However, it is possible to explain the late-time cosmic 
acceleration for the modified gravity scenario in which 
the Lagrangian density is given by C = R + f(G), where 
f(G) is an arbitrary function in terms of Q 1221 . pro- 
vided the function / satisfies some conditions [231 ]. This 
is equivalent to a theory with a scalar field coupled to 
the GB term in the absence of a kinetic term [13, 0, HII . 
A number of f(Q) models that have a matter era fol- 
lowed by a de Sitter (dS) attractor have been proposed 
in Ref. [23| (see also Refs. [2614281 ]). These models can be 
also consistent with local gravity constraints for a wide 
range of parameter space [29( . 

In order to test the cosmological viability of f(G) dark 
energy models, it is important to study cosmological 
perturbations responsible for structure formation. In 
Ref. [27j the evolution of density perturbations has been 
discussed for non-relativistic matter with an equation of 
state w = 0, under the approximation that the back- 
ground cosmological evolution mimics that of the ACDM 
model. In this paper we derive the equation for density 
perturbations with a general constant equation of state 
w. Therefore our analysis includes the perturbations in 
radiation (w = 1/3) as well as those in non-relativistic 
matter. Moreover we use concrete f(G) models that sat- 
isfy both local gravity constraints and cosmological con- 
straints at the background level. This is particularly im- 
portant when we discuss the evolution of perturbations at 
late times, because the deviation from the ACDM model 
can be significant. 

We will show that, in the Universe dominated by a 
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single perfect fluid with an equation of state parame- 
ter w, the perturbations in the fluid exhibit violent in- 
stabilities for w > —1/2 in the small-scale limit. This 
is associated with a negative speed squared c 2 of one 
eigenvector-mode. The perturbations of non-rclativistic 
matter as well as radiation, at some scale, will start to 
grow exponentially during the matter/radiation domina- 
tion, unless the deviation from GR is very small. In GR, 
this same mode does not exist, so that there is no smooth 
limit from one theory to the other. 

This paper is organized as follows. In Sec. [TT] we re- 
view f(G) models that satisfy cosmological constraints at 
the background level as well as local gravity constraints. 
Sec. Hill is devoted to the analysis of cosmological per- 
turbations in f{G) gravity in the presence of a perfect 
fluid with the equation of state w. In this section we 
discuss the presence of an instability at small scales due 
to a negative speed squared for the propagating mode. 
In Sec. IIVI we analyze more in detail for perturbations 
in non-relativistic matter in order to describe the growth 
of large-scale structure. We shall numerically integrate 
the perturbation equations for concrete f(G) models and 
estimate how much deviation from GR can be allowed by 
the observations of galaxy clustering in the linear regime. 
We conclude in Sec. [V] 



II. DARK ENERGY MODELS BASED ON /(g) 
GRAVITY 



Here w e s is an effective equation of state defined by 



Let us start with the following action 
I 



S = 



16ttG 



d 4 xV^U{R + f(g)] + s m , (2) 



where G is a bare gravitational constant and gu is the de- 
terminant for the space-time metric g^. For the matter 
action S m we shall consider a perfect fluid whose equa- 
tion of state w — Pm/ Pm is strictly constant, where p m 
and p m are the pressure and the energy density respec- 
tively. Taking the variation of the action @ with respect 
to gf^j,, we obtain the field equation 



G ui> 



RpyfJ(7(.L RpcrQvfi Rfii^dcrp 



+R fl ag vp + (R/2)(g^g ap - g^g vp )]V p V a f,g 

HQf.g - f)g^ = STrGav , (3) 

where f.g = df/dG, = R^ v - (l/2)Rg f _ lL , is the 

Einstein-tensor, and T^ v is the energy momentum ten- 
sor of matter. 

For the flat Friedmann-Lcmaitrc-Robertson- Walker 
(FLRW) background with a scale factor a, we obtain the 
following dynamical equation 



3ff 2 = Gf,g -I- 24H 3 f. g + 8nG Pr , 



(4) 



where H = a/ a, a dot represents a time derivative in 
terms of cosmic time t, and the GB term is given by 



Q = 2AH 2 (H 2 +H) = -12ff 4 (l + 3w e ff) ■ 



(5) 



Wcff 



-1 - 



2H 
3H 2 



(6) 



The matter energy density p m satisfies the standard con- 
tinuity equation, 



Pn 



3H(1 + w)p rt 



0. 



(7) 



which has the solution p m cx a~ 3 ^ 1+w "> for constant w. 

It is possible to realize a late-time cosmic acceleration 
by the existence of a de Sitter (dS) point that satisfies 
the condition 3H 2 = Gifg(Gi) — f(Gi), where Hi and 
Gi are the Hubble parameter and the GB term at the dS 
point respectively. The condition, 



0<Hff^gg(H 1 ) < 1/384 ; 



(8) 



is required from the stability of the dS point [2J|. We 
have g < and g > during both radiation and matter 
domination. However the GB term changes its sign from 
negative to positive during the transition from the matter 
era {Q = -12H 4 ) to the dS epoch {Q = 24ff 4 ). For 
the existence of standard radiation and matter eras we 
require that f : gg = d 2 f/dg 2 > for g < C?i (23|. Since 
the term 24H 3 f ig in Eq. (g]) is of the order of -ff 8 /,GS, 
this is suppressed relative to 3H 2 for H 6 f^gg <C 1 during 
the radiation and matter domination. In order for this 
condition to hold, we require that f.gg approaches +0 in 
the limit \g\ —J- oo. Recall that even around the de Sitter 
point the condition H 6 f, gg <C 1 is satisfied from Eq. (jSJ). 

A couple of representative models that can satisfy these 
conditions are 12311 



(A) f{9) 



(B) f(G) 



g ( g 

A — == arctan — 

vu: \g* 



1 i— ( G 2 

2 X ^H 1 + g! 



—a\y G 

G ( G 

A arctan — 

VG* \G* 



- aAv G* , 



(9) 
(10) 



where a, A and G* are positive constants. The second 
derivatives of / in terms of G for the models (A) and (B) 

are f,gg = A/[£* 3/2 (l + G 2 /G 2 )] and f,gg = 2A/[£, 3/2 (l + 
G 2 /G 2 ) 2 }, respectively (both of which arc positive for all 

Q)- 

The quantity defined by 



(11) 



is constant for the ACDM model, f(Q) = -2A + cQ (here 
we have included the linear term cG because this also 
gives rise to the equations of motion same as those in the 
ACDM model). In order to discuss cosmological pertur- 
bations in the next section, it is convenient to introduce 
the following quantity 

H = Hi = HGf, gg 

= 72H 6 f,gg [(1 + w eS ){l + 3w otf ) - w' cff /2] , (12) 
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where a prime represents a derivative with respect to 
TV = In a. This quantity characterizes the deviation from 
the ACDM model. During the radiation and matter dom- 
ination one has p = l92H 6 f y gg and p = 72H 6 f y gg, re- 
spectively, whereas at the de Sitter attractor p, = 0. 
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FIG. 1: Evolution of fi (multiplied by 10 4 ) and w e g versus 
the redshift z = ao /a — 1 for the model ([9]) with parameters 
a = 100 and A = 3 x 10~ 4 . The initial conditions are chosen to 
be x = -1.499985, y = 20, and Q m = 0.99999 (see Appendix 
B for the definition of x and y). 

In Fig. [1] we plot the evolution of p and w e g in the 
model (A) for a = 100 and A = 3 x 10~ 4 . In this case 
the quantity p is much smaller than unity in the deep 
matter era (iu ff — 0) and it grows to the order of 10 -4 
prior to the accelerated epoch. This is followed by the 
decrease of p toward with small oscillations, as the 
solution approaches the de Sitter attractor with w e g = 
— 1. For smaller a and larger A, it is also possible to 
realize larger maximum values of p such as p m &x ^ 0.1. 
The qualitative behavior shown in Fig. [T] is generic for 
viable f(Q) models at the background level. 



III. COSMOLOGICAL PERTURBATIONS 

In order to study cosmological perturbations in f{Q) 
gravity we introduce a perturbed metric with 4 scalar 
perturbations a, /?, <j) an d 7 about a spatially flat FLRW 
cosmological background, [3fJ 



ds 2 



-{l + 2a)dt 2 - 2afii dt dx % 

-a(tf [(1 + 2(f>)5 ij + 2d t d 3 j] dxMx j . (13) 



Let us decompose the energy-momentum tensor into 
background and perturbed parts, i.e. T® = — (pm + Sp m ) 



and T° = —p m v ya , where v is a velocity potential. We 
define the gauge-invariant matter density perturbation 

S m , as 



_ Sp m p m 

m = 1 V . 



(14) 



Pm Pn 

We also introduce two gauge-invariant combinations 

$i = 4> + Hv, (15) 
$ 2 = 6£ + iv, (16) 

where <5£ is the perturbation of the quantity £ = fg. 



A. Perturbation equations 

Following a similar procedure to the one developed re- 
cently in Ref. [3l[, one can show that the dynamics of 
cosmological perturbations in f(G) gravity in the pres- 
ence of a perfect fluid with an equation of state w reduces 
to that of two propagating fields $1 and $2 defined by 
the following perturbed action 

5S = J d 4 xL4i<i> 2 + 2^ 2 *i$2 + A 3 <i>l - ,gi(V$i) 2 

-2g 2 V$i • V$ 2 - 33(V$ 2 ) 2 + £($2$1 ~ *1*2) 

(17) 



-m 3 $2 - 2m 2 $i$ 2 ] , 



where At, gi, B, mi are time-dependent coefficients whose 
explicit forms are given in Appendix A. We also have the 
following relation 



a + v = 



1 + ip 
H(l + 6p) 

W c 

-— 0m 

1 + US 



$1 



2H 
1 + 6// 



2H 2 
l + 6p 



(18) 



where p is defined in Eq. ([12)) . 

From the action (fTT|) we obtain the perturbation equa- 
tions in Fourier space 



dt 



-2m 2 $2 - B<& 2 = ; 



dt 



(A 2 $i + ^3*2) + 52^ 2 $i + g 3 k 2 <^ 2 



(19) 
(20) 



+m 3 $ 2 + B$i = 0, 
where k is a comoving wavenumbcr. 

B. Instability at small scales 



Let us show the presence of a small-scale instability 
in f(G) gravity associated with a negative speed squared 
of one propagating mode. This instability appears at 
large redshifts for the models that look like GR (p <C 1) 
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at early times. For sufficiently small scales it is easy to 
show that the highest derivative terms prevail over any 
other terms in the differential equations. Then one can 
approximately find a harmonic oscillator-like dispersion 
relation in the follow way. 

For large k, one can look for an approximate solution 
in the following way. The dominant contribution to the 
equation of motion will be in the form 



.9 V 2 i 



0, 



(21) 



where A is the 2x2 symmetric matrix whose diagonal 
elements are A\ and A3 and non-diagonal element is A^. 
Along the same lines one defines the matrix g with the 
elements g±, g 2 , and 53. Introducing the time-dependence 
<£>j oc exp(iwi) with j = 1,2 in Fourier space, it follows 
that 



(-uj 2 A + k 2 g)<£> w 0, 



(22) 



which has solutions for some values of u> 2 . This approx- 
imation tends to be more accurate for larger k. These 
expressions give the dispersion relation for the propagat- 
ing modes under consideration. 

Non-zero solutions for $ exist provided that the fol- 
lowing relation holds 



det(w 2 A - k 2 g) = . 



(23) 



After finding the eigenvalues w 2 , one can proceed to look 
for the eigenvectors which diagonalize the kinetic opera- 
tor. We find that one eigenvector-mode propagates with 
a speed squared c 2 = w, as expected, and the other one 
with a speed squared 



1 



2H 1 + w 8nGp 7l 
H 2 + l + 4/i 3H' 2 



(24) 



This coincides the result of the vacuum case found in 
Ref. [HI by taking the limit p m — > 0. 

In the Universe dominated by a single fluid one has 
3H 2 ~ 8nGp m and H/H 2 ~ -3(1 + w)/2. Under the 
condition that p <C 1, the speed squared (|24|) reduces to 



-1 - 2w. 



(25) 



This shows the existence of a negative instability for 
w > —1/2. Hence the perturbations in radiation and 
non-relativistic matter are subject to this instability dur- 
ing the radiation and matter domination, respectively. 
In the matter-dominated epoch (f2 m ~ 1) the result (|2"5)) 
agrees with the value c 2 = 1 + 4H/(3H 2 ) that appears 
as a coefficient of the term k 2 /a 2 in Eq. (47) of Ref. [27j |. 
During the transition from the matter era to the acceler- 
ated epoch c 2 , can be quite different from c 2 , because Q m 
is smaller than 1 and the quantity p is not necessarily 
negligible relative to 1. Therefore the stability at late 
times must be checked against the quantity c 2 . 

The reason why the values c 2 and c 2 , arc different can 
be understood as follows. Looking at Eqs. (47) and (48) 



in Ref. [27| . the modes are not diagonalized, as the k 2 /a 2 - 
tcrm still appears for two different fields. Even if the 
coefficient of the k 2 /a 2 term for one of two equations is 
negative, this is not enough to state that there is an insta- 
bility in the system. In other words, even if g± is negative, 
the eigenvalues of A~ x g can still be both positive. An- 
other difference among these two studies is that for the 
evolution of perturbations the authors in Ref. [27[ chose 



a background solution that mimics the evolution in the 
ACDM model 1 . Strictly speaking, this is not a solution 
of the Einstein equations. This can affect the evolution of 
the perturbations especially at late times. In our work we 
use concrete f(Q) models to find cosmological evolution 
of both the background and the perturbations. 

The Laplacian instability mentioned above appears at 
any time in the past at some sufficiently small scales 
(apart from the epoch of inflation during which c 2 , ~ 1). 
This can place tight constraints on f(G) models. Let us 
discuss this more precisely, without any approximation, 
for the growth of large-scale structure during the matter 
domination. 



IV. GROWTH OF MATTER PERTURBATIONS 

Let us focus on non-relativistic matter with the equa- 
tion of state w = 0. In order to derive the equation for 
matter perturbations, we first combine Eqs. (fl~8|). (|19[) . 
([2H and finally take the limit w -+ 0. From Eq. (fTgj) one 
can express $1 in terms of $2, $2 and S m . Multiplying 
Eq. (fT9|) by g 2 and Eq. (j20j) by g\ and subtracting one 
to the other, we find an equation which depends on $1 
only through its first time derivative. Combining these 
two equations and taking the limit w — >• 0, one reaches 
the following equation of motion 



$2 - d 4 <$> 2 



c 2 ^ ) 



where 
cl 
di 
d 2 



1 + il m + 2x + 4p(l + 2x) 
l + Ap ' 
p[l + tt m + 2x + 4p(l + 2x)} 
H(l + An) ' 
fi m [l + 3f2 m + 4x + 4p(l + ix)] 



4(1 + 4^) 



(27) 
(28) 
(29) 



1 Observations show that the current equation of state of dark en- 
ergy is close to —1 to a pretty high redshift, so that the deviation 
from the GR background is not large. Ref. [27ll used this fact to 
neglect the k 2 r) contribution relative to the term already in the 
LHS of Eqs. (47) for the realistic cosmological models one consid- 
ers, for which /i<lor equivalently \H®f gg \ <C 1. In this work 
we make a more general approach of computing the diagonalized 
modes accurately by including the contribution coming from the 
k 2 rj term. 
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d 3 = H 2 {Ax 2 + Ax + x' - 2/i[4(l - 3a; 2 - 2x - x') 

+ 30,„(1 + x) + 8n(2 - 2x 2 - x')]}/[2fx(l + V)], (30) 
3H[1 + O m + 2x + 4/i(l + 2x)] 



dt = — 



1 + 4/i 



(31) 



and x = H/H 2 , x' = x/H and fl m = 8TrGp m /(3H 2 ). 

To find the second dynamical equation we multiply 
Eq. (|T9| by g 3 and Eq. ([20]) by g 2 and then subtract 
the two equations. We then divide it by 5153 — g\ and 
differentiate it with respect to time. This gives rise to 
the equation which involves $1, so that we can replace 
it with 8 m . Furthermore the same equation will contain 
second and third time-derivatives of $2, which can be 
substituted by using Eq. ([26]) . By doing so and taking the 
limit w — > 0, one finds the following dynamical equation 

<L - d 5 8 m - d 6 6 m - d 8 $ 2 + (d 9 + dv^j $2 = , (32) 
where 

2H(l - 2xn + 2aQ , 3H 2 n m (l + x) 

"5 — , , ; , «6 — , , , , (06) 



1 + 4/i 



1 + 4/i 



AH 2 (l + x) , 12H 3 (l + x 

«7 = — — , «8 = - 



(34) 



(35) 



1 + 4/i ' 1 + 4/i 

3H 4 [4:X 2 + 4x + x 1 - 4/i(l - 3x 2 - 2a; - a;')] 
/^(1+4/i) 



Note that in GR /i = and $ 2 = 0. From Eqs. <J30J) 
and (|35| both e?3 and dg diverge in the limit /i — > 0. 
Therefore we can solve Eq. ([26"]) for $2 and substitute it 
into Eq. (|32p . This results in the following equation 

<$ m + dSm + C 2 <5 m = r$ 2 + (d 8 - rrf 4 )$2 , (36) 

where Ci = ?'c?i — ^5, C2 = r^2 — de, and 

M 2 fc 2 k 2 

r=— |, A/l = d 3 + c^— , M 2 =d 9 + d 7 — . (37) 

In order to derive analytic solutions let us consider the 
case fi(k/aH) 2 <C 1. Since we are interested in sub- 
horizon modes (fc 3> aiT), the condition /i <C 1 also fol- 
lows. During the matter domination characterized by 
H ~ 2/(3t) and O m ~ 1 we have 



Ci ~ 2if(l - 2 M ) , C 2 



1 



9 



/' 



fc 

aH 



(38) 

In the GR limit [i ->• and $2 -> 0, Eq. ([351) reduces 
to <5 m + 2HS m — (3/2)H 2 8, n = 0, which has the growing 
mode solution 5 m oc t 2 / 3 in the matter era. 

In the regime ^(k/aH) 2 -C 1 the growth rate of 8 m 
gets larger than that in GR because of the presence of 
the (8/9)/i(fc/aif) 2 term in Eq. ([55]). For sub-horizon 
modes (fc 3> off) this effect is more important than the 
reduction of the friction term C\ induced by 2/i. 



The quantity fi(k/aH) 2 can grow to the order of 1 by 
the present epoch, depending on the wavenumber fc. In 
this case the growth of matter perturbations is signif- 
icantly different from that in GR. During the matter- 
dominated epoch one has d 3 ~ 3H 2 /(2/j.) > and 
c| ~ -1, so that the mass term M\ ~ 3H 2 /(2fi) - k 2 /a 2 
changes its sign from positive to negative at /i(fc/ai?) 2 = 
3/2. This leads to a negative instability for the per- 
turbation $2 through Eq. ([26]) . The evolution of the 
mass term M§ during the matter era is given by ~ 
H 2 (9H 2 / /i— 2k 2 /a 2 ), which changes from positive to neg- 
ative at /i(k/aH) 2 = 9/2. Thus the onset of the negative 
instability can be characterized by the condition 



(aH/ky 



(39) 



In the regime fi(k/aH) 2 ^> 1 one can approximate 
Cjfc 2 /a 2 and M§ ~ d-jk 2 /a 2 , which results in the 
positive mass ratio r ~ 2if 2 /(l + 8/i). Then Eq. ^ 
reduces to 



2H5„ 



H 2 



2(1 + 8/i) 



2H 2 
1 + 8/z 



(40) 



Here we have not used the approximation Notice 
that the coefficient in front of 6 m is positive and hence 
this term does not lead to the growth of 8 m . However, 
the rapid growth of $2 induced by the negative c 2 works 
as a source term for the amplification of 8 m in Eq. (|40[) . 
We have c 2 = 1 and dj = 4H 2 at the late-time dS point, 
which means that both c 2 , and d-? change signs from neg- 
ative to positive during the transition from the matter 
era to the accelerated epoch. Hence we can expect that 
the growth of matter perturbations ends before reaching 
the dS attractor. 

In Fig. [5] the evolution of matter perturbations is plot- 
ted for the model ([9]) with parameters a = 100 and 
A = 3 x 10~ 4 (sec Appendix B for the detail of numeri- 
cal integration). In this case the quantity /i reaches the 
maximum value /i max = 3 x 10 -4 around the rcdshift 
z = 3.6 (see Fig. [1]). Using the criterion ([39]), the per- 
turbation with fc w 60aH is about to enter the negative 
instability region. The quantity [i decreases rapidly af- 
ter it reaches the maximum, whereas aH at z = 3.6 is 
not much different from clqHo today (z = 0). Hence 
one can estimate that the modes with fc < 60ao-ffo are 
hardly affected by the negative instability. In the numer- 
ical simulation of Fig. [5] this can be confirmed for the 
mode fc = 30ao-ffo- Meanwhile Fig. [2] shows that the 
modes with fc > 100ao-£fo exhibit violent negative insta- 
bilities. Note that the apparent discontinuous behavior 
seen in Fig.[2]for the mode fc = 150aoi?o comes from the 
fact that S m temporally becomes negative. 

The wave-numbers relevant to the observed galaxy 
power spectrum in the linear regime corresponds to 
30a o ff o < k < 600a H (i.e. 0.01 h Mpc" 1 < fc < 
0.2hMpc _1 ). For the model parameters used in Fig. [5] 
the resulting matter power spectrum is certainly ruled 
out from the observations of large scale structure. In 
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FIG. 2: Evolution of S m versus the redshift z = ao/a — 1 
for the model (|9]) with the same model parameters as given 
in Fig. [T] We choose three different wave-numbers: (i) k = 
150ao-Ho, (ii) k = lOOao/io, and (iii) k — SOaoHo- The initial 
conditions are x = -1.499985, y = 20, fi m = 0.99999, and 
S m = 5 m /H = W- & , $ = $ = 0. 

Fig. [5] we have chosen the initial conditions $ = $ = 
as a minimal case, but non-zero initial values of $ and 
<E> lead to even larger amplitude of S m . We also note 
that, irrespective of the forms of f(G) models discussed 
in Sec.HU the behavior of perturbations is similar to that 
discussed above. 

The only way to avoid this negative instability is to 
make the parameter fj. as small as possible by changing 
model parameters, so that the modes relevant to the mat- 
ter power spectrum never reach the regime /i(k/aH) 2 = 
0(1). If we take the smallest scale k sa 600ao-£/o of the 
linear matter power spectrum, the condition under which 
the negative instability can be avoided translates into 

/W < 10~ 6 , (41) 

where we have used the approximation aH ss aoHg at 
M = Mmax- Hence the deviation from the ACDM model 
is constrained to be very small. Nonetheless, even by 
introducing by hand this effective cutoff for the wave- 
length due to the experimental apparatus we use to ob- 
serve data, the theory does possess an ultra-violet (UV) 
instability, no matter how small but non-zero fi is. In 
this case perturbation theory at some small scale will 
break down, during anytime in the past up to the dark 
energy domination. Therefore, these theories cannot be 
studied by using perturbation theory, and in general, one 
should expect strong dynamical deviations from GR, as 
the background is not trustable any longer. 

Furthermore we wish to stress that the negative insta- 
bility cannot be avoided as we go to smaller scales. In 



order to avoid violent growth in the non-linear regime 
of the matter power spectrum (k > 600ao-ffo), the con- 
straint on /i max becomes even severer than the one given 
in Eq. f4"Tj) . Moreover the growth rate of matter per- 
turbations gets enormously large for increasing k. The 
point is that we can always find the wavenumber k sat- 
isfying [i(k/aH) 2 w 1 even for very small values of /i. 
This property also persists for the perturbations in ra- 
diation. Since the quantity /i during the radiation era 
is suppressed relative to that during the matter era, the 
scales of instabilities of radiation perturbations are much 
smaller than those of matter perturbations. The only 
way to consistently remove this UV instability is to set /x 
identically equal to zero, that is, the gravitational theory 
exactly reduces to GR. 



V. CONCLUSIONS 

In this paper we have studied cosmological perturba- 
tions in f(G) gravity, in the presence of a perfect fluid 
with a constant equation of state w. In the Universe 
dominated by a single fluid with w > —1/2, we have 
shown the presence of an instability associated with a 
negative speed squared of one eigenvector-mode. Hence 
the perturbations in radiation and non-relativistic matter 
are affected by this instability during the radiation and 
matter domination, respectively. Our results are more 
general than those given in Ref. [27[ in the sense that we 
have considered a general equation of state w and that we 
have not assumed the ACDM-likc background evolution. 

A useful quantity that characterizes the deviation from 
the ACDM model is \i = Hf.g. In the limit that /n0 
(i.e. the ACDM model) one can avoid the appearance 
of the negative instability. If fi ^ 0, the instability of 
perturbations appears for /i > (aH/k) 2 . Even for tiny 
values of \x much smaller than 1, there are small scale 
modes that satisfy this condition. We have studied the 
evolution of non-relativistic matter perturbations numer- 
ically and confirmed that the perturbations are strongly 
amplified once they enter the regime /i > (aH/k) 2 . From 
the requirement that the matter power spectrum in the 
linear regime is not affected by this violent instability, 
we have found that the maximum value of the deviation 
parameter is constrained to be /i max < 1CP 6 . 

Nonetheless the UV limit of this theory remains un- 
satisfactory, as perturbation theory would break down 
eventually at some scale, and the background is not un- 
der control any longer at these scales. When this hap- 
pens, there is not even an easy way, for the cosmological 
background of this theory to be checked against observa- 
tions. It is mostly this UV unpredictable behavior which 
sets the strongest bound. This feature will always re- 
main unless one sets /x identically to at any time, as 
in this case the theory reduces to GR and the unstable 
eigenvector-mode automatically disappears. In this sense 
we believe that this theory is ruled out by our analysis, 
which explicitly shows the existence of eigenmodes with 
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negative squared speed in the past. 

While we have focused on linear perturbations, non- 
linear effects become important once S m grows to the 
order of 1 . It may be of interest to see whether such non- 
linear effects strengthen or weaken the growth of pertur- 
bations. 



Grant-in-Aid for Scientific Research Fund of the JSPS 
Nos. 09314 and 30318802. S.T. also thanks financial 
support for the Grant-in-Aid for Scientific Research on 
Innovative Areas (No. 21111006). 



ACKNOWLEDGEMENTS 



We thank M. Sandstad for discussions and comments. 
The work of A. D and S. T. was supported by the 



Appendix A: Coefficients of the action <| 1 T[) 

Here we present the coefficients appearing in the action (|17|) : 



A l = 

A 2 = 
Az = 

91 = 

92 = 

93 

D 



(12 H 5 £ 3 w + 3 H A wi 2 + 16 H 2 £ 2 tt Gp m w + 16 H 2 £ 2 ir Gp m + 8 tt Gp m H£ + 8 tt Gp m wH£ + wp m n G + tt Gp rn )a 3 

2ttGw(1 + 6£H) 2 H 2 

(-12 H A wi 2 - 3 w£h 3 + 8 tt Gp m Hi + 8 tt Gp m wH£ + 2 it Gp m + 2 wp m TT G)a 3 

2wGtt{1 + 6£,H) 2 ~' 
(12 wiH 3 + 3 wH 2 + 4 wp m TT G + 4 tt Gp m )a 3 H 2 
2wGtt(1 + 6£H) 2 ~' 

(12 H 5 i 3 + 3 £ 2 H 4 + 24 H 3 i 3 H + 24 H 2 i 2 ir Gp m + 6 H 2 H£ 2 + 24 H 2 £ 2 ir Gp m w 



2ttG(1 + 6£H) 2 H 2 

+8 7T Gp m Hl; + 8 TT Gp m wH£ + Wp m TT G + TT G/9 m )d, 

(-12 i 2 H 4 - 3 £H 3 - 24 H 2 Hi 2 - 6 £HH + 2 tt Gp m + 2 wp m ir G)a 

27rG(l + 6££0 2 _ ' 
3aH 2 (4 £H 3 + H 2 + 8 £HH + 2H + A wp m % G + 4 tt Gp m ) 

2ttG(1 + 6^) 2 ' 
(1 + 4 £ff)(-6 w£H 3 + 18 + 3 Hw + 4 wp m 7r G + 4 tt Gp rn )Ha 3 



AwGtt(1 + 6£,H) 2 



m 2 = -B/2 , 



m 3 



4(1 + 6£H) 3 £wGtt 



a 6 H -12 wHH A - 3 wHH + 24 H°w£ + 432 H 'w^ + 192 H°wC 



(42) 
(43) 

(44) 

(45) 
(46) 

(47) 

(48) 
(49) 



192 h£w 2 p 2 m iT 2 G 2 + 72 it Gp m w 2 £H 3 + 240 H 2 £ 2 Hir Gp m + 72 H£Hit Gp m + 384 H£wp 2 m ir 2 G 2 + 144 tt Gp m w 2 £ 2 H 4 
96 tt Gp m w£ 2 H 4 + 80 tt Gp m wiH 3 - 432 H 3 £ 3 H 2 w - 252 H 3 Z, 2 Hw - 288 H 2 £ 2 H 2 w - 48 H 2 ^Hw - 72 H£H 2 w 



864 H 5 Hw£ 3 - 132 H 3 Hw£, - 552 H^Hwt? - 432 H^Hw + 192 H^G'^p 



t3 i, 



TTGp m £ 2 H 4 + 8TTGp m £H* 



+ 672 H 2 i 2 Hwp m TT G + 432 H 2 £ 2 Hw 2 p m iT G + 72 H£_Hw 2 p m Tr G + 144 H£Hwp m TT G - 6 toiP 



(50) 



Appendix B: Numerical integration of dynamical equations 



In addition to the dimensionless variables a; = H/H 2 and fi TO = 8TrGp m / (3H 2 ) we introduce another variable 

y = H/H*, where H* is a constant related to G* (a typical scale of the GB term for dark energy) via H* = QV^. 
Then the background equations can be expressed as 



-4x - ix 



" g/,g ~ / 
H 2 



3(1 -fin,) 



(51) 
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V = X V , 



(52) 
(53) 



The quantities # 6 /,<?e and (Qf t g - f)/H 2 can be expressed by x and y once the model is specified. Introducing the 
following quantity 



the perturbation equations (j2l)]) and (f3"2")l can be written as 



H 

ck 
H 



3x 



-^- + 2 — 
H 2 H 



2^x + 2x 2 - 2x' + c\ 



k 

aH 



$ + diH5' +d 2 S„ 



H 4 H 2 \aH 



2 — -x 

H 3 



4>. 



Numerically we solve these equations together with the bac kground equations (|5ip-(f5"3"l). 
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